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Abstract. We have studied the ^Be(p, 7 )*B reaction in the Halo effective field theory 
(EET) framework. The leading order (LO) results were published in Ref. 1(1] after the 
isospin mirror process, ’Li(n, 7 )*Li, was addressed in Ref. m In both calculations, one 
key step was using the final shallow bound state asymptotic normalization coefficients 
(ANCs) computed by ab initio methods to fix the EFT couplings. Recently we have de¬ 
veloped the next-to-LO (NLO) formalism (3l. which could reproduce other model results 
by no worse than 1% when the ^Be-p energy was between 0 and 0.5 MeV. In our recent 
report p4|, a different approach from that in Ref. (T) was used. We applied Bayesian anal¬ 
ysis to constrain all the NLO-EFT parameters based on measured S -factors, and found 
tight constraints on the S -factor at solar energies. Our S (E = 0 MeV) = 21.3±0.7eVb. 
The uncertainty is half of that previously recommended. In this proceeding, we provide 
extra details of the Bayesian analysis, including the computed EFT parameters’ probabil¬ 
ity distribution functions (PDFs) and how the choice of input data impacts final results. 


1 Introduction 

The ^Be(/7, 7 )®B reaction cross section around zero energy directly affects the solar neutrino flux 0. 
(An expanded discussion on the reaction’s importance for astrophysics can be found in Ref IH.) The 
previously recommended S (0) = 20.8 + 0.7 + 1.4 eV b ©. The first error comes from the capture 
experiment uncertainties. Due to the Coulomb repulsion barrier, the cross section is exponentially 
suppressed, so that the lowest energy of available data is above 0.1 MeV. Extrapolating the data 
down to zero energy based on various models gives rise to the second quoted S (0) error. Therefore 
eliminating the model dependence is critical for improvement. 

Fortunately different scales in this reaction are well separated, and thus the problem is “simple” 
enough to be treated in Halo EFT. The ^Be and proton are considered as two fundamental particles 
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and as a shallow bound state of the two; the high momentum modes resolving the short distance 
(< 3 fm) physics are integrated out, and their effects are subsumed to contact terms in the lagrangian. 
Although the idea is close to the previous potential models’, the EFT puts the scale separation on 
rigorous grounds and should increase in accuracy as terms of higher order in momentum are retained. 
The number of parameters increases while going to higher order, but it is kept minimal, and the error 
of truncation can be estimated. Moreover the important effect of the Coulomb interaction between ^Be 
and proton can be treated exactly. The reaction is entirely an electric-dipole (£1) induced transition 
from the s- and ^f-wave initial states to p-wave configurations in the final ®B. The resulting S -factor 
El is proportional to 


S(E)cxZ,Cl 


|*Sec (E; as, rj + LsSsd (E; Gs, r^) 


e^^Scx {E-,as, rs) f -H \Dec{E)\^ 


The proportionality factor is in Ref. IH. The total spin s of ^Be-p can be 1 or 2, defining two different 
reaction channels. The LO term >Sec is the amplitude of the external direct capture to the ^Be-p 
component in ^B; the same capture but to the lowest excited core ^Be* plus proton component is a 
NLO contribution denoted as esScx, thus the dimensionless factor le^l is small. LsSsy) is another 
NLO term coming from the capture in short distance (< 3 fm); the coefficient Ls ~ 3 fm. The three 
amplitudes in each channel depend on the s-wave scattering length Us and effective range due to the 
initial state interaction effect. Dec is the El transition from the d-wave initial state to the final bound 
state, which is not affected by the initial state interaction at LO. The overall factor in each channel, 
Cj, is the ANC of the two-particle wave function (with total spin s) in ^B. 

There are nine EFT parameters, Cjj, ai, 2 , ^ 1 , 2 , Ei^, and ei (62 = 0 based on angular momen¬ 
tum conservation). Their so-called posterior PDF can be extracted from direct capture data by using 
Bayesian analysis 0. In our analysis, another five parameters ^ 1 , 2 , 3 , 4,5 were introduced to float the 
normalizations of the five capture data sets used. (In total 42 data points with energies up to 0.5 MeV 
were used; Fig. [^provides some information about them.) A Markov chain Monte Carlo (MCMC) 
algorithm with Metropolis sampling 0 was used to compute the 14-dim posterior PDF Q. Note 
prior knowledge about the parameters was incorporated by assigning a prior PDF to each of them. 
For ai^ 2 , it was a Gaussian distribution centered at the measured value 0 with the width equal to the 
experimental uncertainty; the other EFT parameters were given flat prior PDFs over finite value win¬ 
dows 0. The priors were also Gaussian distributions but centered at 0 (the data set’s normalization 
factor is 1 - 1 -^,); the widths were the corresponding experimental common mode errors (CMEs), 2.7%, 
2.3%, 11.25%, 5%, and 2.2% in the order of being mentioned in Fig. (Ref. 0 contains mostly 
relative cross sections, so its CME is handled differently.) 


2 EFT parameters 

Fig. [TJshows a few interesting 1-dim parameter posterior PDFs. (Each PDF is the 14-dim PDF with 
the other 13 parameters marginalized, and was computed via drawing its histogram based on the 
MCMC samples 0 .) We see in the first panel that the data prefer a large C| and a small Cj, and 
the sum is tightly constrained (Cj -H = 0.564(23) fm“* 0 ). The second panel shows the ^Be-p 
s-wave effective ranges. Although we assigned positive flat priors for them, we also required that no 
s-wave resonance exists below 0.6 MeV 0 . The computed PDF of ri disfavors large ri value, while 
that of r 2 ’s is almost flat. (The scattering lengths’ PDFs are close to their priors.) In the third panel, 
each Ls distribution is broad but does have a significant peak around a few fm, overlapping with the 
Cl ,2 ~ 3 fm region suggested by the EFT power counting. Thus we believe that we have detected 
the direct capture at short distance from the data, although the large distance capture dominates the 
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Figure 1. In the first panel, the red solid, blue dashed, and black dashed lines are the posterior PDFs of Cj, C\, 
and C[ + C 2 ; in the second (third) panel, the red solid and blue dashed lines are the PDFs of (Ls) in i = 1 and 2 
channels, respectively; in the last panel, the red dashed (for the Junghans BEl data set), red dotted (the Junghans 
BE3), blue dashed (Filippone), blue dotted (Hammache 1998), and black solid (Baby) lines are the PDFs of the 
corresponding data normalization variable f,. The data sets are shown in Fig.[^ 




Figure 2. The left panel shows the S-factor median values (solid blue curve) and the 1-cr intervals (shaded in 
green) at different energies. The data used in our analysis and a few others above 0.5 MeV are shown, including 
Junghans el.al, BEl and BE3 CSl (filled black circle and filled grey circle), Filippone et.al, im (open circle). 
Baby et.al., I12II13I (filled purple diamond), and Hammache et.al, Ed (filled red box). In the right panel, we 
show three different sets of the S (E) median values and l-tr bounds, as computed by using different combinations 
of the data sets. The blue dashed lines use the combination of Filippone, Hammache and Baby sets, the black 
dashed lines add the Junghans BE3 set, and the red solid lines use all the sets. 


reaction. The last panel shows the PDFs of The corresponding data sets are plotted in Fig.[^ which 
shows minor tension between the recent Junghans data sets and the rest. This is also reflected in the 
PDFs in the last panel; reducing the Junghans BEl and BE3 normalizations is favored, the Eilippone 
and Hammache normalizations need to be increased a bit, and Baby’s is not modified significantly. 
However all floating norms (i.e. 1 + f,), while 1, are entirely consistent with the quoted CMEs. 

3 S factor 

The left panel in Eig.[^ shows the computed 1-cr lower and upper bounds and the median value for 
the S factor. The three quantities were inferred from the S -factor distributions at individual energies, 
which were computed via making the histograms of S based on the MCMC samples. The data are 
plotted for comparison. At the shown energies, the S factor uncertainties « 3.5% or less, and are 
correlated. A numerical table can be found in Ref. a. It is interesting to ask how the choice of data 
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sets impacts our results; this is illustrated in the right panel. Without having the BEl and BE3 sets 
in the analysis, the S factor uncertainty ss 9% at zero energy; adding either BEl or BE3 narrows the 
error bands significantly within the bounds provided by other data. Note that all the three l-cr intervals 
overlap with each other. We also see that combining BEl with the other four sets does not change the 
constraints below 0.1 MeV much; at higher energy the BEl’s impact shows up by constraining S (E) 
further and systematically favoring larger S (E). 

We also repeated our analysis using different priors for the scattering lengths and introducing a 
N^LO term into the capture amplitude, and we found negligible differences. We also found that the 
current 0.7% uncertainty of *B’s proton separation energy could change S (0) by about 0.75%. 

4 Summary 

In summary, the ^Be(p, 7 )*B reaction has been studied in Halo EET up to NLO. We applied Bayesian 
analysis to constrain the EET parameters. The analysis gave a stringent constraint on the total ANC 
squared, and showed that positive short-distance-contribution parameters Li 2 were preferred. Minor 
tensions among the modern capture data sets were seen. We found that S(E) is constrained to ss 3.5% 
between 0 and 0.5 MeV. The S (0) results were robust against exclusion of either Junghans data set 
and against addition of an N^LO term to the reaction amplitude. 
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